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In order to prove the existence of a critical end point (CEP) in the QCD phase diagram it is 
sufficient to demonstrate that at zero temperature E = 0 a first order phase transition exists as a 
function of the baryochemical potential ji, since it is established knowledge from ab-initio lattice 
QCD simulations that at jj. =0 the transition on the temperature axis is a crossover. 

We present the argument that the observation of a gap in the mass-radius relationship for compact 
stars which proves the existence of a so-called third family (aka "mass twins") will imply that the 
T = 0 equation of state of compact star matter exhibits a strong first order transition with a latent 
heat that satisfies Ae/Ec > 0.6. 

Since such a strong first order transition under compact star conditions will remain first order 
when going to symmetric matter, the observation of a disconnected branch (third family) of com¬ 
pact stars in the mass-radius diagram proves the existence of a CEP in QCD. Eor the equation of 
state of the twins the quark matter description is based on a QCD-motivated chiral approach with 
higher-order quark interactions in the Dirac scalar and vector coupling channels. Eor hadronic 
matter we select a relativistic mean-field equation of state with density-dependent couplings. 
Since the nucleons are treated in the quasi-particle framework, an excluded volume has been 
included for the nuclear equation of state at super-saturation density that takes into account the 
finite size of the nucleons. 

Eurthermore we show results of a Bayesian analysis (BA) using disjunct M-R constraints for 
extracting probability measures for cold, dense matter equations of state. This study reveals that 
measuring the radii of neutron star twins has the potential to support the existence of a first order 
phase transition for compact star matter. 
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1. Introduction 


Neutron stars are compact objects which serve as astrophysical laboratories for studying ex¬ 
tremely dense matter in their interiors |[T], possibly composed of exotic forms of matter such 
as hypernuclear matter and deconfined quark matter [Q]. Their investigation can be performed from 
both sides: astrophysical observations and nuclear collision experiments. Observational programs 
aim at measuring properties like mass, radius, moment of inertia, spinning and cooling rates among 
others. Nuclear experiments in the other hand look for precise determination of properties of nuclei 
in the laboratory via nuclear collisions and reactions. In the last ten years important developments 
have occurred which allowed to place constraints to the equation of state (EoS) of neutron stars 
iQ]. These measurements are, however, not sufficient yet to completely determine the structure and 
composition of these compact stars and many questions remain open. The recent observation of 
the most massive neutron stars of about 2 M 0 ^ has resulted in both ruling out EoS models that 

cannot support such massive stars and in restricting the range of densities in compact star matter. 

Theoretical studies of dense matter have explored the QCD phase diagram and have shown 
that a smooth crossover transition occurs at very low chemical potentials and finite temperature - 
conditions for which colliders like the EHC and RHIC have provided compatible data. At higher 
chemical potentials and sufficiently low temperatures a first order phase transition to deconfined 
matfer mighf occur which implies fhe exisfence of af leasf one critical end point (CEP) separating 
both phase transition types. The most massive neutron stars are likely to harbour deconfined quark 
maffer in fheir cores [Q ^ where values of fhe baryon density well exceed fhaf of fhe nuclear 
safurafion densify and isospin asymmetry is extreme with proton fractions below 10%. The ex¬ 
istence of a CEP and onset of deconfinement (CPOD) results in important features in the neutron 
properties under certain conditions. The so called mass twin phenomenon being a most pronounced 


consequence of this [10, 11, 12]. It is therefore relevant to point out what those observables are and 
to assess the feasibility of their detection. In this contribution we present the recent theoretical de¬ 
velopments in the study of massive neutron stars related to the high-mass twin phenomenon which 
can provide evidence for the existence of the CEP in the QCD phase diagram. 


2. The neutron star twins 


A strong first order phase transition to deconfined maffer in neufron sfar matfer will resulf 


in a disconnected branch in fhe mass-radius relafion [ ]13| ] often called fhe fhird family afler white 
dwarfs, and neufron sfars have been infroduced as firsf and second families, respecfively. The 
compacf objecfs on fhe fhird family branch have a similar mass buf smaller radii fhan pure neufron 
sfars wifh a difference of up fo fwo kilometers for fhe models discussed here. This properly led to 
the name "mass twins". While being known in the literature for a long time [ [T^ , [Ts] , p^ , |T7| , the 
possibility of their existence at the high mass of 2 M 0 was demonstrated only recently ©[ni’B]- 
A new microscopic approach for the description of hybrid stars was recently introduced in [ [T^ ] 
incorporates important ingredients for the neutron star EoS and can also provide the mass twins. By 
introducing the effect of excluded volume of baryons due to the Pauli exclusion principle between 
the quarks contained inside each nucleon the hadronic EoS stiffens up. On the other hand multi¬ 
quark interactions in the quark matter description also stiffen the EoS at high densities [p^. 
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This results in a hybrid EoS for dense quark-baryon matter which supports massive hybrid 
stars of more than 2M0, eventually as a third family of compact stars. That case is equivalent to 
the existence of high-mass twins, with considerable radius differences. 


2.1 Why do we care about the twins? 


The suggestion to prove the presence of a CEP in the QCD phase diagram by the possible 
existence in nature of high-mass neutron star twins is in principle testable by astrophysical ob¬ 
servations. Besides this striking aspect of high-mass twins, their existence could solve also the 


presently most urgent problems of compact star structure |20]: 


The masquerade problem. The EoS of hybrid stars with a smooth phase transition can look 
almost the same as the pure hadronic one. This results in a degenerate mass-radius relation 


which masquerades hybrid stars as hadronic ones []21[]. 


The hyperon puzzle. The appearance of hyperons in neutron star interiors softens the EoS 
limiting the maximum stable star mass configuration that in most of the cases is much below 


the detected 2M0 value [22]. 


The reconflnement puzzle. In certain EoS a deconfinement transition occurs at lower den¬ 


sities but at much higher values the hadronic EoS becomes favourable again [123|,|24{]. 


The measured radius issue. The radius determination by measurements of thermal spectra 


from quiescent low-mass X-ray binaries inside globular clusters in [ 25 j is one of the lowest 


values reported. On the contrary, both measurements of hot spots in | |2q ] and phase-resolved 
spectroscopy studies report rather large radius values. The nature of massive neutron 
stars twins feature favor long radii providing a resolution to the radius issue. 


The neutron star twins features are rather robust against variations of the EoS, e.g., due to con¬ 
sideration of structures in the coexistence phase of the first-order phase transition ("pasta phases"). 
Although the appearance of pasta phases in the hadron-quark interface will have an effect on the 
mass-radius diagram it will remove the third branch only for strong effects as explored phenomeno¬ 
logically in 1^]. The third branch will persist also for fastly rotating stars as shown in where 
stability against radial oscillations was investigated. Studies for the presently discussed high-mass 
twin hybrid star EoS show that fast conversion scenarios from neutron stars to hybrid stars, e.g., 
by mass accretion are possible also under fast rotation for objects like the high-mass millisecond 
pulsar PSR J1614-2230 [^] with a period of 3.15 ms. 


3. Bayesian Analysis for hybrid neutron stars 


Bayesian analysis is a powerful statistical tool that can provide probabilities of detection of 
the hybrid twins phenomenon. In the context of constraining the EoS with observational data for 
masses and radii of compact stars it has first been introduced by [p0|, 31]. However their studies 


using X-ray bursters suffer from uncertainties related to the determination of the stellar atmospheres 
which is a crucial ingredient in spectral studies. Moreover, the parabolic region in the M-R diagram 
which follows from the determination of the luminosity radius is largely degenerate with an 
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EoS of the APR type and thus fakes a high probability for such an EoS. In [32, 331 the authors 
take a different approach by using a flexible parameterization of the EoS capable of including twin 
configurations on the one hand and is based on different neutron star measurements on the other, 
not relying on the atmosphere composition and being disjunct in mass and radius determination. 


4. Future studies on the twins 


There is still space left for exploration concerning the neutron star twins. Erom the stellar 
evolution side the open question of the feasibility of a dynamical transition from a hadronic neutron 
star into the hybrid one via accretion. As mentioned before, rotation of hybrid stars has already 
been explored in therefore predictions for the highly rotating massive neutron stars are in order. 
The physical imprints of gravitational radiation from excited modes is also a matter of investigation. 
Previous studies for hybrid stars going into this directions can be found in [3^ 35]. 

A future precise, simultaneous measurement of both a short and a large radius value for two 
objects is a crucial element for twin identification. Eictitious inclusion of such a measurement into 
the Bayesian analysis studies shall provide more realistic probabilities for the twin phenomenon 
appearance. 


5. Conclusions 


In this contribution we have reviewed the mass twin star phenomenon as a possible solution 
to several problems in modern astrophysics: the masquerade problem, the hyperon and reconfine¬ 
ment puzzles and the neutron star radius issue. Moreover, the observation of high-mass twins in 
nature would prove the existence of a critical point in the QCD phase diagram that is currently 
being looked for in terrestrial experiments with heavy ion collision experiments. Since simulations 
of lattice QCD at finite temperature and vanishing baryon densities find a crossover behaviour, the 
observation of high-mass twin stars implying a strong first order phase transition at zero temper¬ 
ature and finite baryon density would be strong support for the existence of at least one critical 
endpoint in the QCD phase diagram. We have also pointed out that Bayesian analysis is a tool ca¬ 
pable of assessing the feasibility of high-mass twin detection. In concluding we express our hope 
that future observational programmes with satellite missions like NICER [|^] or terrestrial large 


scale observatories like SKA [ )37[ ] will provide the necessary accurate simultaneous determinations 
of masses and radii for a sufficiently large set of neutron stars which possibly will find support for 
the high-mass twin star phenomenon and thus the QCD critical endpoint. 
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